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Abstract 


Currently, there is a need for a high-speed, high-lift civilian transport. Although 
unconventional, a flying wing could fly at speeds in excess of Mach 2 and still retain the capacity 
of a 747. The design of the flying wing is inherently unstable since it lacks a fuselage and a 
horizontal tail. The project goal was to design, construct, fly, and test a remote-piloted scale 
model flying wing. The project was completed as part of the NASA/USRA Advanced 
Aeronautics Design Program. 

These unique restrictions required us to implement several fundamental design changes 
from last year's Elang configuration including wing sweepback and wingtip endplates. Unique 
features such as a single ducted fan engine, composite structural materials, and an electrostatic 
stability system were incorporated. The result is the Banshee '94. Our efforts will aid future 
projects in design and construction techniques so that a viable flying wing can become an integral 
part of the aviation industry. 
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1. Introduction 


Currently, there is a need for a high-speed, high-lift civilian transport. Although 
unconventional, an oblique flying wing could fly at speeds in excess of Mach 2 and retain the 
capacity of a 747. The long term goal of this project is to produce a scale model of such an 
aircraft. This project, sponsored by NASA/USRA, dealt with the design, construction, flight, and 
testing of a remote-piloted scale model conventional flying wing. It is hoped that design and 
construction techniques learned from this project will aid in the future construction of a scale 
model remote-piloted oblique flying wing. 

1.1 The Flying Wing 

The flying wing was first developed in the 1940’s by Northrop Aircraft Inc. . During the 
World War II, military planners wanted a bomber that would be able to fly from the East Coast 
of the US to Germany, round trip, with payload. This goal was defined as "10000 pounds for 
10000 miles” 6 . Jack Northrop believed that a flying wing would be able to attain this goal. He 
believed that the lack of a fuselage and tail would reduce weight and, therefore, produce a higher 
cruise lift to drag ratio. Also, the lower weight would make it possible for the wing to carry 
more fuel. These design considerations would give the flying wing the range necessary for an 
intercontinental bomber. The two bombers constructed, the B-35 propeller driven and the B-49 
jet propelled flying wings (see Figure 1.1), were inherently unstable due to the lack of a tail. 
Aerodynamic forces produced unfavorable pitching moments in the flying wing. To correct this, 
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the flying wing is swept back with its control surfaces behind the center ot gravity, producing 
some stability. A damping system was also added to account for the pitching moment of the 
wing. 



Figure 1.1 The Northrop B-35 and B-49 Flying Wings 16 

In the end the Convair B-36 and not the Northrop B-49 that was chosen as the next 
intercontinental bomber. In November 1949, the existing B-35’s and B-49's were destroyed. 
Derivatives of the Flying Wing design can be seen today in the Northrop B-2 Stealth Bomber 
(see Figure 1.2) and the recently cancelled McDonnel Douglas A-12. 



Figure 1.2 The Northrop B-2 Stealth Bomber 10 



1.2 The Oblique Flying Wing 


An oblique flying wing differs from a conventional wing, like the B-35 and B-49, since 
it, flies with one w'ing tip ahead of the other, relative to the direction of flight. The concept of 
the oblique wing was first developed after World War II by NASA aerodynamicist, Robert T. 
Jones. This design, due to the orientation of the wing, would reduce drag increasing the lift to 
drag ratio at any speed. The high lift to drag ratio of the oblique wing aircraft would improve 
low speed performance for an aircraft designed to tly at high speeds 1 '. 

The first oblique wing demonstrator aircraft was the NASA AD-1 (see Figure 1 .3). This 
design placed a high aspect ratio wing on a pivot atop a conventional fuselage. The wing itself 
could be swept up to 60 degrees. The AD-1 was first built in 1978 but work was halted soon 
after due to the fuel crisis. The 50 test flights of the AD-1 proved highly sucessful and showed 
that such technology could be adapted for a flying wing configuration 1 * 



Figure 1.3 The Oblique Wing demonstator, the NASA AD-l y 

Airlines and aircraft manufacturers are showing serious intrest in the oblique living wing 



as a high speed, high capacity commercial transport. The conceptual oblique flying wing (see 
Figure 1.4) could carry 500+ passengers, more than a conventional Boeing 747, and fly at speeds 
in excess of Mach 2. It is believed that an oblique flying wing would be 16-30% cheaper to fly 
than the Boeing 747. Travelling twice the speed of a 747, the oblique flying wing would use 
slightly more fuel while carrying more passengers. 

The obliqueness of the wing would be changed by the rudders. Orientation of the wing 
would determine the speed of the aircraft; the lower the angle, the lower the speed. Overland 
flights could be made just below the speed of sound with high efficiency just by changing the 
angle of the wing. 

The eventual goal of these NASA/USRA sponsored projects in WPI's Advanced 
Aeronautics Design Program is to construct a scale model oblique flying wing. For the last two 
years, the scope included designing and constructing a conventional scale model flying wing. 
The results are the Elang and the Banshee '94. The current, the Banshee '94, design improved 
upon and corrected the mistakes made in the design of the Elang. Also, the Banshee '94 used 
base numbers of the Elang for calculation and scale of the aircraft. 

1.3 The Flying Wing, Elang 

Last year's design resulted in the Elang (see Figure 1.5). The Elang used the custom 
HPRS 33 airfoil (see Figure 1.6). To ensure minimum induced drag (elliptical distribution), the 
Elang used an elliptical planform. The flying wing had a root chord of 45.5cm, a span of 2.5m, 
and an aspect ratio of 7. To provide yaw control, two SD 8020 symmetrical airfoil (Figure 1 .6) 
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were used as vertical stabilizers. 



The propulsion system for the lilang consisted of two ducted fan engines. Ducted fans 
were chosen over propeller driven engines because of the ducted fan's clean aerodynamics. 
Propellers tend to produce undesirable turbulence over the wing and control surfaces, adding to 
the instability of the aircraft Each of these ducted fan engines operated at 22000rpm and 
produces 15.6N of thrust. Due to their small size and high complexity, the ducted fan engines 
were reluctantly replaced with standard propeller driven engines late in last year's design. 
Propeller engines are highly accepted in the remote control aircraft world and are easier to 


operate and maintain 



2) Jirfoi I vec'ion 



PA. 11 ; airfoi I wr[,on 


Figure 1.6 The HPRS 33 and SD 8020 airfoils 

The materials used in the construction of the Elang needed to be strong enough to 
withstand "accidents” and have a low weight. The Elang used a styrofoam core with a skin of 
carbon fiber. The carbon fiber skin proved to be strong enough for the testing stage of the 
Elang . 

The testing stage for the Elang was very brief. The Elang's lack of an automatic stability 
system proved to be fatal to the aircraft. The Elang was fitted with a standard R JC gyroscope 
system, but it was unclear whether it was used during test flights. Without some form of 
automatic stability system, a flying wing, on any scale, will not fly. 

1.4 The Flying Wing, Banshee f 94 

The purpose of this project was to design, construct, fly, and test a scale model remote 
piloted flying wing. As in industry, an emphasis was put on design teams. The project group 
was split into four separate sub-groups, each dealing in a different aspect of design, working 
together towards the final goal. The result, the Banshee 1 94 (see Figure 1.7), named after the 
spirit of Celtic lore, was designed to solve the problems involved with the stability of flying 
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wings. The Aerodynamics group dealt with calculation of lift and drag, also with airfoil 
selection. The Controls and Stability sub-group was responsible for design and analysis of 
control surfaces for the Banshee '94 and also the installation and construction of an automatic 
stability system which will control the inherently unstable design. The Propulsion group dealt 
with the selection and testing of an engine with suitable power for the flying wing. The 
Structures and Construction group dealt with material selection, finite element modeling of that 
material, and the actual construction of the flying wing. 



Figure 1.7 The Flying Wing, Banshee '94 





2. Design Summary 


2.1 Introduction 

Before the detailed calculations and methods are discussed, a basic summary of the 
aircraft’s design will be given. In this section, the general design parameters for the Banshee 
'94 will be discussed. A table comparing the design parameters of the 1992-1993 Flying Wing, 
the Elang, to the Banshee ’94 is given. 

2.2 General Design Specifications 

The Banshee '94 is a tapered, swept flying wing. This differs from the Elang , in that 
the Elang was an elliptical flying wing. The Banshee '94 was designed with a straight wing for 
ease of construction and a tapered, swept wing for increased stability. A four view computer 
aided design (CAD) drawing of the Banshee ’94 is shown in Figure 2. 1 . The aircraft is powered 
by a single ducted fan engine located in the center of the wing. The Elang was designed with 
two ducted fan engines, which promoted stability problems due to the difficulty in equalizing the 
two engine’s thrust. Each wing tip on the Banshee ’94 has an endplate to act as a vertical 
stabilizer. For control surfaces, each wing has an aileron and a elevator and each endplate has 
a rudder. These will enable control in all three axes; roll, pitch, and yaw. 

2.3 Specific Characteristics 

2.3.1 The Wing 

The total wingspan is 92.0 inches. The root chord is 18.0 inches and the tip chord is 8.0 



Figure 2.1 The Banshee '94 4-view CAD drawing 


inches. The leading edge is swept back 15 degrees from the horizontal. This provides an aspect 
ratio of 6.689, a taper ratio of 0.444, and a total planform area, including the duct, of 
approximately 8.5 square feet. From the analysis, the lift to drag ratio for the Banshee ’94 is 
19.44. 


2.3.2 The Propulsion System 

The Banshee ’94 is powered by a 4.6 horse-powered OS Max 91 scaled-aircraft engine. 
The engine has a thrust of 11.0 to 14.0 pounds-force from manufacturer’s data and a thrust to 
weight ratio of 4.67. The engine powers a nine-bladed ducted fan unit which propels the 
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aircraft. The engine is fueled by two fuel tanks enclosed inside the wing on either side of the 
duct. 


2.3.3 Electronic Components 

The engine and control surfaces in the Banshee ’94 are controlled by electronic servos. 
These servos are controlled by wireless remote. To insure stability, an electro-static stability 
system was built and installed inside the aircraft. All of the components are housed within the 
wing and accessible through removable panels. 


2.4 The Cruise Condition 

This report details the analyses performed which led to the final design of the Banshee 
’94. These analyses resulted in cruise conditions for the aircraft. The aircraft will fly with the 
) wing at an angle of attack of four degrees to insure stability. The aircraft will have a cruise 

velocity of thrity-three miles per hour to insure enough lift to remain in level flight. 


2.5 The Banshee ’ 94 Versus The Elang 


Wingspan 
Root chord length 
Wing sweep 
Aspect Ratio 
Total planform area 


The Banshee [ 94 
92.5 inches 
18 inches 

15 degrees 


The Elang 

98.4 inches 
18 inches 
0 degrees 


6.689 

8.5 square feet 


7 

9.6 square feet 


) 


1 1 



Lift to drag ratio 


19.44 


19.69 


tt of engines 

Total thrust from engines 
Angle of attack 
Cniise velocity 


1 

14 Ibs-f 
4 degrees 
33 miles per hour 


2 

7 lbs-f 
4 degrees 
34 miles per hour 


) 


) 
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3. Aerodynamics 


3.1 Introduction 

Each individual design group is responsible for designing and analyzing specific aspects 
of the Banshee ’94. The Aerodynamics group is responsible for determining the aircraft’s lift 
and drag after determining the aircraft’s configuration and airfoils. In addition to building and 
calibrating the electro-static stability system, the Controls/Stability group must determine the 
stable angle of attack and effect of control surface deflections. The Propulsion group chooses 
the engine and analyzes its thrust at various flight and throttle conditions. They must also 
calculate take-off and landing parameters including take-off speed and take-off distance. The 
Structures group is responsible for structural analysis, material selection, and fabrication of the 
aircraft. 

Although each group possesses significantly different responsibilities, each group must 
perform its design and analysis in cooperation with the other groups. To obtain cruise 
conditions, the Aerodynamics, Propulsion, and Stability groups need to perform a combined 
analysis. With the chosen airfoils, the correct statically stable angle of attack may be determined 
by the Stability group. Once the angle of attack is determined, the Aerodynamics group 
determines the velocity, lift, and drag at cruise conditions. With the drag known, the Propulsion 
group can determine the correct engine throttle setting to produce the necessary thrust at cruise 
conditions. 

For structural analysis, all groups needed to provide information to the Structures group. 
The Controls and Propulsion group provided component weights while the Aerodynamics group 



provided the lift and drag distribution. This information allowed for the use of a finite element 
analysis of the flying wing. 

Also, every group needed to work in cooperation while determining the lofting and 
construction. The Aerodynamics group needed to a configuration that would prove easy to 
manufacture and an airfoil with the proper thickness for component placement. The Stability 
and Controls group needed to provide control surface and component sizes to allow the 
Structures group to construct the wing. The Propulsion group needed to provide engine size so 
that a properly dimensioned duct could be manufactured. With this general overview, further 
discussion of the design process is discussed below. 


3.2 Aircraft Configuration 



Figure 3.1 The Banshee 94 with quarter chordlines 


The configuration of the Banshee ‘94 incorporates a twisted swept wing with endplates. 
Wing sweep creates two major effects: decreased drag and additional stability. Decreased drag 
improves an aircraft’s performance. The stability effects are created by the airfoil’s 
aerodynamic centers acting along the length of the root chord axis rather than a single point 
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(quarter chord), as is the case with an unswept wing (see Figure 3.2). The specific stability 
calculations will be discussed in another section. 

The Banshee ’94 incorporates aerodynamic wing twist in its design. Twist varies the 
sectional lift coefficient along the wingspan allowing for increased stability. Two types of wing 
twist can be utilized: geometric, where the wing itself is physically twisted to change the 

sectional angle of attack: and aerodynamic, where the airfoil section changes along the wing. 
The Banshee '94 flying wing uses the latter. 

Another feature of the Banshee’s configuration are endplates at the end of both wing tips. 
Being the largest wing tip device, this particular wing tip device is rarely used due to their 
excess skin friction drag produced by the considerably large surface area. However, on a flying 
wing, endplates may be used as vertical stabilizers. On conventional aircraft, (with tail and 
fuselage) the vertical stabilizer is usually placed on the tail. On such an aircraft, wing endplates 
would only serve as a lift to drag improvement device. By allowing the endplates to serve a 
dual purpose, little additional drag is created. This placement allows the best performance of 
the vertical stabilizer. On conventional aircraft, the vertical stabilizer is placed on the tail. On 
an aircraft without a tail or fuselage, the tips of a swept wing are the rear-most locations on the 
aircraft. 

In addition to acting as a control surface, the endplates also increase the lift to drag ratio. 
The endplates allow for a non-zero lift distribution at the wing tip (see Figure 3.3). In addition 
to increased lift, the endplates also reduce drag. Large plates placed at the end of the tips of 
the aircraft inhibit the trailing vortices which create induced drag. The aerodynamic effects 
influenced by the endplates increases overall performance. 



Figure 3.3 Lift distribution on wing with and without an endplate 
3.3 Airfoil Selection 

Due to the use of aerodynamic twist, a number of airfoils are necessary. The Banshee 
’94 uses a series of five airfoils, the Eppler 325 series 5 . The E325, E326, E327, E328, and 
E329 airfoils are pictured in Figure 3.4. The airfoil sections are placed on the wing as 
illustrated in Figure 3.5. The aerodynamic data for the E325, E327, and E329 is listed in 
Appendix Al. 
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Figure 3.5 Airfoil placement along the 
wing 


For a flying wing, certain considerations require additional attention during design, such 
as stability and component placement. The chosen airfoils must have certain characteristics 
including necessary moments, thickness, and similarity of shape. The E325 series satisfies these 
requirements. 

With the absence of a tail, a horizontal stabilizer cannot be incorporated in the design. 
Without such a device, an aircraft could not achieve pitch stability with a conventional airfoil. 
Most airfoils have a negative moment about the aerodynamic center. With the horizontal 
stabilizer absent, it becomes necessary to use a set of airfoils which produces a positive pitching 
moment about the wing. To obtain a positive moment, an airfoil with a reflexed camber is 
needed. For a swept wing, the reflex of the camber increases along the wingspan toward the 
tip, as illustrated in Figure 3.4. This change in the reflex of the camber creates aerodynamic 
twist along the wing promoting pitch stability. The airfoils near the root, the E329 and E328, 
have significant negative pitching moments. However, with proper placement of the airfoils, 
the lift produced by the airfoils creates a positive moment about the aircraft’s center of gravity. 
The calculation of center of gravity is discussed in depth in the controls and stability section. 
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In addition to the absence of a tail, a flying wing has no fuselage. This poses an 
additional problem with component placement since all components must now be placed within 
the wing. There must be adequate space for the components. The E325 series has one of the 
largest thicknesses of all available airfoil series for tailless aircraft. 

A series of airfoils implies a similarity in shape. In fact, the mid- and semi-span airfoils 
are geometrically linear extrapolations of the tip and root airfoils. This similarity allows for 
additional benefits during analysis and construction. The aerodynamic data varies linearly 
among airfoils, as illustrated in Figure 3.6. This linear relationship allows for certain 
assumptions to be made about averaging the airfoils aerodynamic characteristics. This will be 
discussed in detail in the aerodynamic and stability analysis sections. Another benefit occurs 
with construction. The hotwiring technique used allows the form of each half of the wing to be 
cut swiftly and accurately. The similarity of shape allowed the wing to be fabricated using only 
two templates at the root and tip. This construction method will be discussed further in the 
structures section. 

3.4 Drag Calculations 

3.4.1 Introduction 

The main purpose of this section is to describe in detail the calculations performed to 
predict the drag on the Banshee ’94. The first section describes the different assumptions used 
for the analysis. The second section gives a detailed description of the calculations; citing all 
of the required equations and constants. 
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Figure 3.6 C, vs. C m from Airfoil Design and Data, Richard Eppler 


3.4.2 Assumptions 

As with any calculations, certain assumptions needed to be made in order to analytically 
determine the drag exerted on the Banshee ’94. These assumptions are summarized as follows: 

1 . Incompressible Flow: Since the Banshee ’94 will fly at a cruise velocity well 
below the speed of sound, compressibility effects are negligible. 

2. Turbulent Flow: The historical trends of the cruise velocity of similar aircraft 
demonstrate a low velocity profile. However, the characteristic length of the 
Banshee ’94 is large enough to raise the value of the Reynold’s number to about 
700,000. This value is well above the 500,000 transition Reynold’s number. This 
is important to the skin friction calculations. 


) 
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3. 


> 


Flat-Plate Theory: All of the major components of the aircraft including the 

landing gear and endplates were treated as flat-plates with the same cross-sectional 
area. The landing gear were assumed to be one inch thick and three inches tall. 
The wing was not treated in this matter. Experimental data from the airfoil 
reference provided the drag coefficients. 

4. The Wing Sizes: The wing was "divided" into different sections representing the 
different airfoils. The chord length was held constant at an average of the chord 
lengths associated with each section. 

5. Constant Properties of Air: Since the Banshee '94 will be flying at low altitudes, 
the properties of air remain constant with slight changes in altitude. Sea 
level values were assumed. 

6. A pproximate Pressure Drag: Since the only way to accurately account for 

pressure drag is through a lengthy Computational Fluid Dynamics procedure, an 
approximation was used as a more efficient method. This was only implemented 
in the calculations for the endplates. 


3.4.3 Calculations 

The following constants were used in the analysis: 


density of air 

P 

0.076474 lbm/ft 3 

viscosity coefficient of air 

P 

1.2024 x 10 5 lbm/s ft 

Oswald Efficiency Factor 

e 

0.98 


) 
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A component method was used to calculate the drag exerted on the Banshee '94. The 
first component that will be considered is the actual wing itself excluding the landing gear, duct, 
and endplates. The wing was sectioned analytically to represent the three airfoils that data was 
available for; E325, E327, and E329 (see Figure 3.5). The drag was calculated for each section 
and then summed. 

By definition, the drag is calculated as follows: 

Drag=C D j pV 2 S 

For the wing, the parasitic drag was calculated for each portion of the wing using the different 
drag coefficients for each airfoil (see Appendix Al). These values were then summed to provide 
a total parasitic drag for each wing. The induced drag is the drag due to the lift. It is calculated 
using the following equations: 

pV 2 S 


C„ r kC 2 L 

k=— 

nAe 

As with the parasitic drag, each induced drag component was calculated for each section of the 
wing and then summed for a total induced drag for each wing. 

Thin-plate theory was used to calculate the drag on each of the endplates. The coefficient 
of drag for one endplate is calculated using the following equations: 
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where S WET is the surface area of the endplate and S REF is the reference area (planform area) of 
the endplate. Using the previous equations, the drag was calculated for each endplate. 

The drag for the duct was calculated in a similar manner. The coefficient of drag for 
the duct was calculated as follows, neglecting pressure effects: 
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where is the surface area of the duct, Sref is the reference area (planform area) of the duct, 
and C F is calculated the same way as for the endplate. 

The drag for the landing gear was calculated in the same manner. The coefficient of drag 
for the landing gear was calculated as follows, neglecting pressure effects: 


C D =c d =c n 

U LC U o w U F ba 


‘ WET 


U Flg 5 

^ REF 

where Swet is the surface area of one landing gear, S REF is the reference area (planform area) 
of one landing gear, and C F is calculated the same way as for the endplate and the duct. 

The total drag is calculated by adding up all of the drag from all of the components: 
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Drag to =2 *Drag w +2 * Drag E + Drag D +3 *Drag LG 


W - wing; E - endplate; D - duct; LG - landing gear 

3.5 Lift Calculations 

3.5.1 Introduction 

The lift for the Banshee ’94 was an important consideration. Using an estimated weight 
of 13 pounds, it was necessary for the cruise conditions to generate 13 pounds of lift to 
overcome the weight. 

3.5.2 Calculations 

The following section describes the calculations used to calculate the lift at different 
cruise velocities. The only constants used in these calculations are the density of the air, the 
different coefficients of lift obtained from published experimental data, and the different 
reference areas for the sections of the wings. Since it was assumed that the cruise altitude would 
not be significant enough to warrant a change in the density, a standard sea level density was 
used; 0.076474 lbm/in 3 . The coefficients of lift and the reference areas can be found in the 
spreadsheet in Appendix A2. 

Each wing was divided into three sections to represent the three different airfoils for 
which coefficient of lift data was available. An average chord length was used for each section 
to determine the planform area needed for the calculations. 

The lift for each section was calculated using the following equation: 
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Lift=C L ^pV 2 S 

Each lift was calculated and then summed to provide the total lift exerted on the Banshee ’94 
by the wing. The lift calculations were iterated over a velocity interval to determine the 
necessary speed to keep the Banshee ’94 airborne. 

3.6 Lift Versus Drag Calculations 

3.6.1 Introduction 

One of the major factors in considering an aircraft’s performance and efficiency is the 
lift to drag ratio. The advantage of a flying wing in this area manifests itself in the high lift to 
drag ratio. This phenomenon is due simply to the high planform area of the lifting surfaces 
versus total planform area of the aircraft. For a flying wing, the ratio of lifting surface area to 
total area is nearly one. Below is a summary of the calculations for the lift to drag ratio. 

3.6.2 Calculations 

The lift and drag were determined from empirical methods which were explained in 
sections 3.3 and 3.4 of this report. The lift to drag ratio was calculated using the results of 
these sections. 

3.7 Results and Conclusions 

Below, in Figure 3.7, is a graph of the calculated lift in pounds versus velocity in miles 
per hour. It can be clearly seen from the graph the approximate cruise velocity necessary to 
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maintain a steady altitude is thirty-three miles per hour. This velocity is reasonable for an 
aircraft of this size. The Elang had a cruise velocity of thirty-four miles per hour. 

The spreadsheet that was used for the calculations can be found in Appendix A2. Using 
cruise condition of angle of attack of 4° and thirty-three miles per hour, the total drag exerted 
on the Banshee ’94 is calculated to be 0.617 pounds-force. This value is reasonable for an 
aircraft of this size and shape. The propulsion system being used should provide sufficient thrust 
to achieve cruise conditions. Figure 3.8 is a graph of the total drag force exerted on the 
Banshee ’94 in pounds versus the cruise velocity in miles per hour. 
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Figure 3.8 Drag vs. Velocity 


From these results, the lift to drag ratio for the Banshee ’ 94 was calculated to be 19.44. 
This value is greater than that of a typical aircraft, as would be expected. However, this value 
is close to those found for other such flying wing aircraft such as the Elang. Below is a table 
of aircraft and their subsequent lift to drag ratios: 


Aircraft 

L/D 

Banshee ’94 

19.44 

Elang 

19.69 

Surya NASA Advanced Aero MQP 91-92 

19.55 
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Concord (Trans- Atlantic) 

7.5 

XB-35 Flying Wing 

20 

Oblique Wing Aircraft 

9.8 

B-47 

17.2 

AVRO Vulcan 

17.0 


Figure 3.9 shows the thrust and drag versus velocity curve at different throttle settings. 
The thrust data is detailed in Chapter Five. The figure shows that at twenty-five percent power, 
the thrust equals the drag at the predicted cruise velocity of thirty-three miles per hour. 



Figure 3.9 Thrust and Drag vs. Velocity 





4. Stability and Controls 


4.1 Introduction 

The primary responsibilities of the stability/controls group was to determine the important 
flight parameters to ensure stable flight. The information provided by controls/ stability is 
particularly important because of the intrinsic instability of the flying wing. The most important 
parameters such as the center of gravity, the flight condition angle of attack, and the static 
margin, are all crucial to the success of the flying wing. 

Stability refers to the ability of a vehicle to return to its original equilibrium position, 
without pilot assistance, after a disturbance has pitched, yawed, or rolled it to a different angle 
or has caused a speed change. 15 Static stability is present if the forces created by the disturbed 
condition displace it in the correct direction to return the aircraft to its original position, while 
dynamic stability pertains to the dynamic behavior of the motion of the aircraft before and after 
a disturbance. Dynamic stability is present if the dynamic motions of the aircraft will eventually 
return it to its initial state. 

The rigid body motions of aircraft may be conveniently divided into two classifications, 
longitudinal and lateral motions. Longitudinal motions occur in the plane of symmetry, which 
remain in its original position. Lateral motions, such as rolling, yawing and sideslipping, 
displace the plane of symmetry. 13 The technical significance of this distinction is that for normal 
symmetrical aircraft with small displacements these two types of motion are independent of each 
other. 
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In performing the stability analysis for Banshee ’94, we considered the three primary axes 
of rotation: pitch, roll, and yaw, however our major concern was stability along the pitch axis. 
In conventional aircraft, a horizontal tail is placed so that it may offset the moments originated 
by the wing and fuselage. Consequently, a flymg wing without a horizontal tail demands that 
a supplementary method for stabilization be used. 

4.2 Static Stability 

The goals of the static stability analysis were to find: 1) A workable cruise condition 
angle of attack for Banshee '94\ and 2) A center of gravity for the wing which ensured static 
stability. Our primary objective was to focus on the stability of the aircraft in the pitch axis, 
which is anticipated to be more problematic than the other two axes. Stability in the yaw and 
roll axes is expected to be acceptable because of the symmetry of the wing in these axes. 

The first objective was to derive the classic static stability curve: C mcg vs a. The starting 
point was the lift, drag, and moment coefficients as determined by the aerodynamics group. 
(Table 4.1) 

The moment coefficient about the center of gravity is arrived at through the equation: 



Next, the lift and drag moment arms, which vary with changes in a, the angle of attack 
had to be considered. (Figure 4.1) 
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Table 4.1, Experimental Airfoil Data 5 




E329 

E327 

E325 

N 

a 

c, 

c d 

c 

'-'mac 

c, 

c d 

c 

mac 


c d 

c 

''aw c 

2 

0.53 

0.0086 

-0.05 

0.35 

0.008 

0.0 

0.15 

0.0077 

0.05 

3 

0.65 

0.0091 

-0.51 

0.48 

0.0085 

-0.0024 

0.27 

0.0079 

0.048 

1 4 

0.73 

0.0094 

-0.52 

0.56 

0.0087 



0.0081 

0.046 

5 

0.86 

0.0102 

-0.53 

0.65 

0.0094 

-0.0072 

0.45 

0.0085 

0.044 

r 6 

0.95 

0.0108 

-0.54 

0.76 

0.01 

-0.0096 

0.59 

0.0094 

0.042 


1.03 

0.0115 

-0.55 

0.87 

0.0107 

-0.012 

0.69 

0.0101 

0.04 
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The final equation used is as follows: 
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The derivation of this equation had to account for three different airfoils of the wing due 
to the fact that each airfoil had a different aerodynamic center (because of sweep-back) as well 
as a different planform area. Here, x is the center of gravity minus the aerodynamic center, 
which was assumed to be the quarter chord, and c & S respectively are the mean chord length 
and planform area of that particular airfoil. 

An iterative process was employed through the use of a spreadsheet in order to determine 
a center of gravity of 7.43 inches behind the front of the duct and an angle of attack of four 
degrees for cruise conditions. The C mcg vs a curve for this center of gravity was derived (Figure 
4.2). 

Static stability is ensured due to the negative slope of the plot. If the airplane is 
momentarily pitched to a higher angle, the pitching moment change will tend to restore the plane 
to its original attitude. If the airplane is pitched nose-down, the resulting moment change will 
tend to raise the nose. 15 

This spreadsheet was also used to determine the neutral point and the corresponding 
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Figure 4.2 Static Stability Curve 


static margin of the wing. To arrive at this, C, ncg was plotted against C^,., the average lift of the 
three airfoils. When the slope of this curve is equal to zero, this is the point on the plane which 
has been labeled the neutral point. The neutral point of Banshee '94 was found to be at 8.7252 
inches behind the front of the duct or 1 .2952 inches behind the center of gravity. The graph can 
be seen in Figure 4.3: 

Note that while the slope of this curve doesn’t appear to be zero, the average slope of 
this curve is zero. By placing the neutral point farther towards the front of the wing, the C m vs 
C, curve slopes upwards while placing the neutral point farther back causes the plot to slope 
downwards. 


) 
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The static margin is defined as: 


Static Margin = 


X-X. 


eg 


c 



Where X n is the neutral point and C is the mean chord length. The static margin for 
Banshee ’94 was found to be 0.0996 or 9.96 percent, an acceptable number for conventional 
aircraft. 


) 
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4.3 Dynamic Stability 


The dynamic stability analysis entailed analyzing the wings reactions to sudden 
disturbances placed upon the wing during flight conditions, particularly along the pitch axis. A 
major concern with Banshee '94 is a low moment of inertia about the pitch axis due to the 
absence of a horizontal tail. A horizontal tail would be used to provide damping effects so that 
the wing could retain an acceptable amount of controllability. 

The governing equations of this analysis are: 

M cs =m 

Where M cs is the moment caused by the deflected control surface, found through the use 
of the inviscid vortex panel code, and 5 is the angle of deflection of the control surface. Also 
used is the equation: 


/ 

Where M lol is the sum of the moments about the center of gravity derived from the wing, 
duct, etc., as well as the moment derived from the deflected control surface. I is the moment 
of inertia about a particular axis, found through a solid modeling program, ARIES, and /? is the 
angular acceleration of the wing along that axis. Extremely high or low values for (3 denote poor 
controllability. By observing the angular accelerations, dynamic stability could be determined. 
Plots of vs 5 were arrived at for the pitch, roll, and yaw axes. (Figures 4.4, 4.5, 4.6) 
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Figure 4.4 Angular Acceleration, Pitch Axis 


- ) 

4.3.1 Results 

®max * n the roll axis was found to be around 70 rad/ sec 2 , indicating an acceptable amount 
of control. It can be seen in the plots for the elevators & rudders that the angular accelerations 
in the pitch and yaw axes are small, denoting low controllability. This isn’t a foreseeable 
problem because of the symmetry in the yaw axis, however the absence of a tail in the pitch axis 
__ demands that further measures be taken. To conclude, it was determined that the wing has poor 

dynamic stability along the pitch axis. 


) 


35 


QHK3JNAL PAC* ft 

or poor fuAUTY 




This analysis helped us to calculate the sizes for the control surfaces. Although our 
control surface size estimation was based on Elang’s values, we decided to slightly oversize the 
elevators, ailerons and rudders to ensure better control. The final sizes appear in Table 4.2: 

The dynamic stability analysis ascertained that several alternative design methods would 
have to be employed so that the wing would be stable during flight. Some of these design 
implementations are: 1) Sweeping the wings back for improved stability in all three axes; and 
2) The use of an automatic electrostatic stability system. 
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Figure 4.6 Angular Acceleration, Yaw Axis 


Table 4.2, Control Surface Sizes: 


Control Surface 

Percent of Chord Length 

Percent of Planform Area 

Elevator 

18.0 

13.0 

Aileron 

18.0 

2.9 

Vertical Stabilizer 

N/A 

12.0 

Rudder 

N/A 

31.7 percent of 
endplate planform 
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4.4 Sweep-back 


One feature we decided to implement for the aircraft was to sweep-back the wings. From 
a stability point of view, sweep-back is important for three reasons: 

1) Swept-back wings have large dihedral stability even when they have no dihedral angle. 
This arises because the effect of the angle of yaw is to increase the sweep-back angle of one 
wing panel and decrease it for the other side of the airplane. The change in sweep alters the 
effective dynamic pressure, increasing the lift on one side, lowering it on the other, and 
producing a restoring rolling moment. 13 

2) Sweep-back increases the moment of inertia along the pitch axis, helping to decrease 
the angular acceleration created by moments along this axis. 

3) Sweep-back likewise aids in yaw stability. Figure 4.7 shows that a yawing moment 
produced by the differences in drag tends to restore the airplane to a straighter flight path. 

Even with sweep-back, a flying wing is inherently unstable along the pitch axis. 
Normally, the horizontal tail would offset the moments generated by the fuselage and the wing. 
Because Banshee ’94 lacks a tail, it was decided that another method for stability would have 
to be utilized. 

4.5 Electrostatic Stability System 

The electrostatic stability system, (ESS) was invented by Maynard L. Hill of the John 
Hopkins Applied Physics Labs. It is a device that operates by measuring the earth’s magnetic 



Sweepback yields improved 
Yaw Stability 


Figure 4.7 Effect of Sweep-back 

field as a vertical reference, since the earth's electromagnetic field gradient is about 650 volts 
per foot. 7 A complete two-axis system, functioning in the roll and pitch axes, weighs about one- 
tenth of a pound, uses 80 milliwatts of power and costs about $130. Its principal component is 
the StaticMaster 1 a one inch square strip of radioactive polonium w'hich measures static voltage. 
The alpha radiation that it emits is quite harmless and has a half-life of about six months. 
StaticMasters are also widely used in industrial and hobby applications for the elimination of 
electrostatic charges. 

In a simplified sense, an electrostatic autopilot is a solid-state electronic unit that can be 

1 Supplied by NRD Corporation, 2937 Alt Blvd. Grand Island. NY. 
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used in place of complex inertial gyroscopes whenever the electric field is vertical, a condition 
that on average exists about 90 percent of the time over 90 percent of the earth’s surface. 7 
Electrical disturbances near thunderstorms cause malfunctions over distances that range out to 
about 10 kilometers from the edge of the storm. Also, the atmospheric electric field is sometimes 
useless over broader regions covered by frontal snow, rain, and sand storms. The method is not 
appropriate for use on manned aircraft where reliability is a prime requisite, and this is why the 
concept has not been seriously pursued for use on military weapons and aircraft. 7 

The StaticMaster displaces charged ions from the earth’s magnetic field, and the signal 
from these ions are sent through a differential voltmeter (Figure 4.8). Next, the signal is sent 
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to mechanical servos which stabilize the aircraft accordingly. A more detailed view of the 
controls layout can be observed in Figure 4.9. 
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Figure 4.9 Electric Component Flowchart 


A. Differential Voltage 
Amplifier 

B. Receiver 

C. Elevator Servo 

D. Aileron Servo 


E. Rudder Servo 

F. Throttle Servo 

G. Switch Servo 

H. Switch 

I. Receiver Battery Pack 


J. ESS Battery Pack 

K. Op. Amp. (Gain = -1) 

L. StaticMaster 

M. Servo Interface Board 


This diagram represents the overall interaction between the servos, receiver, and the 
differential voltmeter, (DVA). A mechanical two channel switch, operated by a servo is 
employed here. When the switch is in its left position, the analog signal from the DVA travels 
first through a servo interface board, which converts the analog signal to a servo-compatible 
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pulse-width signal. The signal then moves the control surfaces through the servos, automatically 
stabilizing the aircraft. When the switch is in its right position, a direct link exists between the 
servos and the receiver, so that the control surfaces can be displaced manually. 

4.6 Component Placement - Planform 

Two factors were kept in mind while the planform analysis was executed, center of 
gravity position and ease of wiring. The components and their appropriate weights can be seen 
in Table 4.3. 


Table 4.3, Component Weights and Sizes 


Component 

Mass, lbm 

Volume, in 3 

Density , lbm/in 3 

ESS Battery Pack 

0.1647 

2.4079 

0.0684 

Receiver Battery 
Pack 

0.2093 

3.0596 

0.0684 

Servos (A) 2 

0.0664 

1.1989 

0.0554 

Servos (B) 3 

0.1106 

1.9968 

0.0554 

Fuel Tank (Full) 

0.6644 

23.5590 

0.0282 

Receiver 

0.0808 

3.5459 

0.0228 

StaticMaster 

0.0083 

0.1709 

0.0486 

Switch 

0.0245 

0.8107 

0.0302 

Circuitboard 

0.1909 

10.9688 

0.0174 


ARIES w as used to ensure that the center of gravity was 7.43 inches behind the front of 
the duct. Figure 4.10 displays the final overall planform of Banshee ’94. 


2 Futaba Model FPS148 

3 Futaba Model FPS3001 
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4.6 Conclusions 


From the outset, the stability analysis was performed under the shadow of intrinsic 
instability. The static stability analysis proved that with the center of gravity placed at 7.43 
inches behind the front of the duct, Banshee ’94 would be statically stable due to the C mct vs a 
slope being negative. The dynamic stability analysis, as expected, showed that the deficiency of 
a horizontal tail would present a major design problem. 

This plight was rectified through two pathways. First, sweeping back the wings w’as a 
simple measure which provided immediate stability in all three axes. Finally, placing an 
electrostatic stability system provided a more advanced system in which Banshee ’94 could be 
stabilized in both the roll and pitch axes. 

As anticipated, our calculations showed that the wing would be inherently unstable along 
the pitch axis, sweep-back included. The electrostatic stability system has proven to be a 
working device and we fully expect it to be operational. Sweep-back along with the ESS will be 
welcome additions to our aircraft that will empower us with better control. 
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5. Propulsion System 


5.1 Objective 

In the process of designing Banshee 94 ’ it was imperative to select a propulsion 
system which is suitable for the unique characteristics of a flying wing. As in all aircraft, 
there must be a light weight engine that will produce enough thrust suitable for runway take- 
offs and stable horizontal flight. In our circumstances the propulsion system should also 
conform to the shape of the wing to reduce drag and give the appearance of a flying wing 
(eliminate a fuselage). One of our desired factors was the ability to have excess power to 
recover from any flight difficulties. 

5.2.1 Engine Selection 

The engine selected for our task proved to meet the requirements expected for level 
flight and maneuverability. The engine selected was a single OS Max 91 VR-DF 



Figure 5.1 OS Max 91 VR-DF Engine 
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gas engine (figure 5.1), used in many high performance remote controlled aircraft. 

Along with the OS Max 91, we selected a ducted fan instead of a propeller. The reason for 
this decision was to eliminate propeller wash and to best conform to the wings shape. The 
engine is positioned 4.5 inches from the center of gravity along the centerline of the aircraft. 

5.2.2 Manufacture’s Data 

The OS Max 91 VR-DF Ducted Fan Gas Engine is a relatively new engine on the 
market. It is used in many remote controlled aircraft, such as F-16 and Sabers, to simulate jet 
engines. It is used for this task because it can be placed within the body of the aircraft and 
also because of the high exhaust velocities achieved by this engine. 

The OS Max 91 engine is a two stroke gas engine made especially for ducted fan use. 
A ducted fan engine is designed differently from conventional propeller engines. Ducted fan 
engines perform at speeds of 22,000 RPM, 12 causing the engine to over heat during shorter 
periods of time than conventional engines. Because of this problem ducted fan engines are 
designed with large heat sinks to cool the engine at high RPM’s. The OS Max 91 engine 
weights 25.77 ounces n without the fan attached, and approximately 3.0 pounds with the 
entire assembly. The engine runs on a standard fuel for two stroke engines with 20% 
nitromethane and 20% lubricant. The high percentage of nitromethane enables the engine to 
run at high RPM’s and enhances high performance abilities. This enhancement of 
performance gives the engine a total power output of 4.6 BHP and a static thrust of 14.0 lbs. 12 
This engine surpasses the net power of the previous system and insures that there will be 
sufficient power to maneuver the aircraft. 
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5.3 Methods of Propulsion 

Two main methods of propulsion were considered for the aircraft, propellers and 
ducted fans. Both methods are currently in use for remotely piloted aircraft, but propellers 
are extensively used while the ducted fans are just beginning to become accepted. 

5.3.1 Propellers 

Free propeller driven aircraft are the most widely used propulsion systems in the 
remotely piloted aircraft industry. Propeller engines are very reliable and take very little 
maintenance. The propeller can be easily mounted on to the driving shaft of the engine with 
a nose cone. The simplicity of the engine makes it useable for many aircraft. 

There is one major disadvantage to the use of propellers on some aircraft. On all 
propeller engines there is a turbulent flow of air which is produced by the high speed rotation 
of the propeller, called propeller wash. Propeller wash originates from the engine and travels 
over the airfoil and past the control surfaces, mainly the elevators. The turbulent air over the 
elevators and ailerons makes the control of the aircraft difficult in many situations. Making 
the controls delayed or ineffective for a period of time. This lack of control is undesirable for 
a flying wing. Since a flying wing is an inherently unstable configuration, the lack of control 
would make the aircraft even more unstable. This undesired effect could possibly create an 
error in the performance of the elevators or ailerons, resulting in catastrophic consequences. 


5.3.2 Ducted Fan 
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Ducted fan engines have only recently been accepted as a form of propulsion for 
remotely piloted aircraft. Ducted fan engines are primarily used in scaled jet aircraft with 
internal engines. The reason for their application in these aircraft is because of their ability to 
be mounted inside the frame of the aircraft. 



Figure 5.2 Ramtech Ducted fan Assembly 

Ducted fan engines are high performance engines which run at high speeds and high 
temperatures. The high performance of the ducted fans makes them have a greater output 
than other two stroke propeller driven engine in its scale. One drawback is the fact that 
ducted fan engines are very sensitive due to their high performance, and have been known to 
overheat when not operated with the proper equipment. 

Besides the possibility of overheating, the advantages over propeller driven engines 
pointed our team to work with a ducted fan. The Ducted Fan assembly being used with the 
OS Max 91 engine is made specifically for the OS Max 91. The fan unit is 5.25 inches in 
diameter with nine blades which substitute for a propeller. The nine fan blades rotate 
counter-clockwise and drive the flow of air into ten stator blades which directs the air to flow 
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parallel to the duct. The fan unit is mounted at the neck of the engine to the stator blades, 
and integrates itself to the engine to make an aerodynamic shape for parallel flow over the 
engines heat sinks. The assembly is mounted to the front of a 20 inch duct which passes 
through the center line of the aircraft. The entire assembly has a total weight of 
approximately 3.0 pounds and produces about 14 pounds of thrust. Giving more than 
adequate thrust for the purposes of our aircraft. 

5.4 Single Engine 

Last years flying wing design, The Elang, used two ducted fan engines. The engines 
used last year were the OS 25 Ducted Fan Gas Engines. These engines were placed 22 inches 
from the center line of the aircraft and 1 1 inches from the quarter chord. Each engine was 
predicted to produce 15.6 N (3.5 lbs) of thrust at static conditions. With the ducted fan 
attached, each engine assembly weighed approximately one pound and operated at a maximum 
speed of 22,000 RPM. Both engines received fuel from a central fuel tank located near the 
center of gravity of the aircraft. The central fuel tank was a 12 oz fuel tank, giving the 
approximate flight time of 10 minutes. The twin engine system was difficult to operate. In 
many cases, it was very difficult to run the engines at equal thrust levels. These conditions 
resulted in yaw and control problems. 

As a result the current propulsion team decided to utilize a higher horsepower single 
engine. A single engine will eliminate the possibility of these problems reoccurring. It will 
also be simpler to adjust to the engine during preparations and flight tests. The possibility of 
unequal thrust will be eliminated with only one power source. The fuel mixture can be easily 
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adjusted by the needle valve to regulate the maximum power and fuel consumption. 

5.5 Location 

One of the major factors in the design is the location of the propulsion assembly. The 
engine is the heaviest object on the aircraft, and has the largest effect on the center of gravity 
of the aircraft. After determining the correct location of the center of gravity, the internal 
control components and fuel tanks can be placed accordingly. Then the engine can be place 
to counteract the moments of the internal components. After placement of the internal 
components the engine assembly was placed 4.5 inches away from the center of gravity. This 
location is only several inches in front of the leading edge of the airfoil, and makes the duct 
of the airplane 20 inches in length. 

5.6 Fuel System 

The fuel system is composed of two 12 ounce tanks located on both sides of the duct. 
Twelve ounces of fuel will give the Banshee 94' an estimated engine running time of 20 
minutes. This will give the aircraft adequate time to start the engine, taxi, and fly, without 
concerning about the amount of fuel remaining. The fuel tanks were placed on opposite sides 
of the center of gravity of the airplane. The positioning of the tanks is critical because they 
are the second heaviest internal component, and because of their unstable center of gravity. If 
the fuel tanks are not positioned correctly there is the possibility of reducing the stability of 
the aircraft. The fuel tanks are connected to the engine by a series of tubing which passes 
through the duct and attaches to a Y-connection on the left side of the duct, and then proceeds 
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to the engine inlet. This configuration enables the engine to draw from both tanks 
simultaneously. Both tanks will not need to be pressurized from an exhaust tube because the 
engin e already produces an appropriate vacuum to draw the fuel. 

5.7 Engine Type 

The current engine was selected after a series of experiments using similar engines of 
lower power outputs. The engines were tested in the WPI wind tunnel to determine the 
effects of varying wind velocities on engine output. 

In order to record the data from the experiments, several input devices were utilized to 
determine the variables, thrust, wind velocity, torque, and RPM. The constants of the tests 
were temperature, pressure, and fuel mixture. The input devices consisted of a strain gages, 
the wind tunnel pressure gages, and photocell tachometer or a fiber optic sensor. The engines 
were mounted on an aluminum mount. The mount was made of an aluminum rod which was 
attached to the engine and secured by a stand. The rod was attached to a flexible bar 
containing two pairs of strain gages to measure the thrust and torque of the engine. Each pair 
of the strain gages were connected to voltage meters which amplified the signal in to useful 
data. The thrust gage was calibrated before the experiments. A linear relationship between 
voltage and thrust was determind. The wind tunnel pressure gage was used to determine the 
ffeestream velocity during the experiment. The photocell tachometer and fiber optic sensor 
were used to measure engine RPM’s during engine testing to determine when the engine is 
running at maximum power. 

The tests were conducted on two engines available to the propulsion team. The 
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engines were, the OS 61 propeller engine and the OS 25 ducted fan. Each engine, running at 
max power, was tested at wind tunnel speeds from 0 to 100 mph. A resulting Thrust vs. 
Freestream Velocity plot was determined for each engine. The OS 61 engine had a maximum 
thrust of 5.5 lbs, which is slightly lower than the manufactures data of 6.0 lbs. The OS 25 
ducted fan engine test resulted in 2.0 lbs of thrust, which was lower than the manufactures 
data of 3.5 lbs. The resulting curves of each engine proved to be very similar in shape. The 
experimental results were then scaled up to the current engines performance level using a 
maximum thrust of 14 lbs given by the manufactures data for the OS Max 91 VR-DF ducted 
fan engine. 

In many cases the manufactures data is greater than the actual output of the engine. 

The output of the engine can fluxuate due to an improper mixture of air to fuel ratio or 
improper fuel choice. In the case of a improper air to fuel ratio, it takes practice to realize 
when the engine is working at its maximum RPM. Even though these mistakes may be 
corrected, the manufactures data may still not be accurate. For this reason the propulsion 
team calculated the thrust output of the engine at several engine speeds, 100% (22,000 RPM), 
75% (16,500 RPM), 50% (1 1,000 RPM), and 25% (5500 RPM). The major assumption made 
for these calculations is that the engines have similar properties and react similarly to the 
change in velocity. 

The first calculation was taken at one quarter (5,500 RPM) the power output of the 
engine. The purpose for this calculation is to verify the cruise velocity of the flying wing. 
During flight the engine will not be running a full throttle, but at about 25% to 30% its 
maximum output. This will ensure a reasonable cruise velocity which corresponds to the 
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figures calculated by the aerodynamics team. 

The second series of experiments determined the thrust at half throttle (1 1,000 RPM). 
The main purpose of this calculation was to see the engine output at normal flying conditions, 
because a large percentage of the flight time will be within these limits. 

The third experiment was taken at 75% output (16,500 RPM). This was to see the 
engines abilities during climbs and aerial maneuvers. The fourth calculation was taken at full 
throttle. These experiments to made to see the performance of the engine at the manufactures 
maximum thrust data. 

Plans for tests on the current engine are now in progress. These tests were delayed 
because the current stand was not sufficient to handle the OS Max 91 engine’s increased 
weight and power. The experiment will be run in a similar manner as the smaller engine tests 
were conducted. The proposed tests will coincide with the calculations previously taken, by 
performing the tests at 100%, 75%, 50% and 25% power. 

The results from the engine testing and calculations can be seen in the Thrust - Drag 
vs. Velocity plot in figure 5.3. 

At the intersection of the drag curve and the 25% thrust curve is the cruise velocity of the 
aircraft. The cruise velocity of the Banshee 94 ’, through the engine testing, has been 
estimated at 35 mph. The maximum velocity, assuming the manufactures data is accurate, is 
estimated at 93 mph. The maximum velocity at 75% power during horizontal flight is 
estimated at 76 mph. 
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Thrust * Drag Vs. Velocity 



Figure 5.3 Thrust - Drag vs. Velocity 
5.8 Excess Power 

After reviewing the results of the current engines predicted performance, it is easy to 
notice the excess power delivered by the engine. One of the factors considered when 
choosing the OS Max 91 ducted fan engine was the ability to use extra thrust when needed. 
The extra thrust could be used for a variety of reasons. 

When taking off from a static position, it is necessary to have a greater velocity than 
required for cruise conditions. Therefore cruise conditions can not be the maximum thrust of 
the engine, but only a small percentage of the available power. Take-off velocities are not 
significantly larger than cruise velocities, but the excess power is required in several other 
areas. 

An additional need for power is required in the rate of climb of the aircraft. If a 
engine is selected to perform at only cruise velocity, the aircraft will not be able to climb. 
The greater the engines power the faster and higher an aircraft can climb. When selecting an 
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engine for the Banshee 94 \ a main concern was the flight performance. We wished to give 
the pilot less limitations on the planes capabilities, and more room for error. This was 
another reason for extra power. If the pilot was to make an error in judgement, he would 
have greater chance to compensate for his mistake if there is more power available. The 
additional power would also allow the pilot to perform a larger array of aerial maneuvers. 

Excess power was not only considered for flight conditions, but to compensate for 
errors in calculations. To compensate for the possibility of errors the designer must use a 
material or power source which is greater than the calculations limitations. Some of the 
possible errors in the Banshee 94 ’ may have occurred in the assumptions made for the drag of 
the flying wing. There is also a large possibility of an error in the manufactures data on the 
maximum power output of the OS Max 91 engine. The engine must also compensate for 
unpredicted weight not estimated in the original design. 

One additional factor not noticed in the early stages of design were the carbon fiber 
rod, used to connect the wings and as spars, and the control rods passing through the duct. 
These objects pass through the duct and partially block the air flow produced by the fan unit. 
This problem is solved by creating a symmetric airfoil shell which will contain the control 
rods and carbon fiber rods. The shell will be placed in the duct over the areas of interest to 
reduce the blockage and to make an easy transition for the airflow. 

5.9 Rate of climb 

As discussed previously, the rate of climb of the aircraft was one of the major reasons 
for utilizing an engine with excess power. From the test data the propulsion 
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group was able to make calculations on the rate of climb at various power outputs and 
velocities. These equation used to predict the rate of climb of the aircraft was 


RC = Climb rate T = Thrust V = Velocity W= Weight D = Drag 


RC = 


60 * 


( TV-DV} 

W 


The results formed a plot of climb rate vs. velocity for power outputs of 75% and 50%, as 


shown in figure 5.4. 



Figure 5.4 Climb Rate vs. Velocity 


The maximum rate of climb with 75% thrust is 19 ft/s at 55 mph. 


5.10 Equations 


Other calculations which resulted from the Thrust - Drag vs. Velocity plot were, lift- 
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off velocity, lift-off distance, and landing distance. The lift-off velocity estimates the speed at 
which the aircraft is able to take off with a 4 degree angle of attack and without using the 
elevator control surfaces. 

p = density W= Weight C, = Coeficent of lift 
S = Surface Area V s = stall speed 

Lift-off Velocity 

K, = 1.2V, = 1.2 

= 56. 2 ft Is = 38. 81 mph 

As discussed be for the lift-off velocity is slightly greater than the cruise velocity of 33 mph. 

The lift-off distance was also determined using the same assumptions of a 4 degree 
angle of attack and without the use of elevator surfaces. 


2W 

pSC, 


T = Thrust D = Drag g = gravity 
Lift-off Distance 


D b 


1. 44 


pSC, 


2 W 
2( T-D> 
Wg 


Thrust at 75% 


D h = 64 . 59 ft 


The estimated landing distance is predicted using a velocity at fifty feet of 40 mph and 

a landing velocity of 35 mph. 

V L = Landing Velocity a = acceleration 

d - A A 

D dect! 2 a la D *ce! = 59 . 83 ft 
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5.11 


Conclusion 


After considering the needs of the Banshee 94 a light weight, high power, single 
engine with a adequate rate of climb, stable flight, and reasonable running time. The 
propulsion team decided upon the OS Max 91 DF-VR Ducted Fan Engine. The OS Max 91 
engine satisfies the needs of the Banshee 94 and leaves room flight or design changes. The 
OS Max 91 engine was also purchased to provide for next years Advanced Aerodynamics 
Design project, expecting that they will not increase the design parameters significantly. The 
engine meets the expectations of the current Advanced Aerodynamics Design group and will 
hopefully continue to meet the expectations of future design groups. 
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6. Structures 


6.1. Introduction: 

The structures group, working with the Aerodynamic group, the Control group 
and the Propulsion group surmounts the problems of maintaining the structural integrity of 
Banshee '94 The structures group is responsible for the structural integrity of the aircraft 
The structural integrity has to be accomplished with the lightest possible weight since any 
excess weight has a detrimental effect upon the aircraft performance Another 
responsibility of the group is to determine the geometric shape of the structures, the 
mechanical loads, the deformations acting through the structural boundary, and the 
material make-up of the structure The deformations have also to be determined to assure 
that there is no excessive displacement existing in the structure, since a large deflection 
may make the plane stall and become unstable The structures group analyses on stress 
not only to determine the kind of material to use, the appropriate thickness, the size and 
the cross sectional shape of the plane, but also to assist in designing of connections for all 
parts of the aircraft Ultimately, the group is to determine the best structural arrangement 
for items such as power plants, built in fuel tanks, and landing gear, and to respond to the 
important requirement of the structural design safety with light weight 

In order to accomplish these tasks . as well as to obtain a detailed analysis for a 
given aircraft structure, the structures group uses two different Finite Element pakages I- 
Deas and MSC/Nastran The structures group used simple analytical methods for 
confirming the results from the Nastran Program; especially with regard to the 
displacement (or deformation) of the wing In addition, the structures group conducted 
experimental tests on test beam section made of PVC foam covered on both top and 
bottom with Kevlar This section was tested to failure; the data obtained was compared 
against the values given by the FEM codes 
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6.2. Material selection: 

For construction materials, there are two basic categories; conventional 
construction, and composite construction Conventional construction (balsa wood, 
aluminum etc ) can lead to time consuming fabrication when applied to a complex wing 
shape They are also brittle in nature, and damage easily if not handled with care 
Composite construction, however, can afford great strength, and needs little attention for 
proper maintenance Comparing these two factors, we decided that a composite structure 
of a face sheet, or "skin", laminated to a foam "core" would provide the best strength to 
weight ratio as well as the desired durability for the aircraft 

In the early 1960's, the emergence of boron filaments started a new generation of 
composites Composites with modulus continuous filaments like boron and graphite are 
referred to as advanced composites This remarkable class of material has had a 
profound impact on airframe design technology An advanced composite material is 
defined as a material consisting of small-diameter, high strength, high modulusfor 
stiffness) fibers, which is embedded in an essentially homogeneous matrix. This makes the 
material anisotropic, since such materials have mechanical and physical properties that 
vary with direction In addition to advanced composite materials such as Kevlar, and new 
matrix materials including polymide, thermoplastics, and even metals such as aluminum, 
titanum, and nickel are also used Using the properties of composite materials, the weight 
is reduced by 30%. Even though the resulting structures are more expensive than their 
metal counterparts, the composite materials are still ideal for the structural applications 
because of their important properties, high strength-to- weight and stiffness-to-weight 
ratios Among the new materials, Kevlar has been chosen for the "Banshee 94". 
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6.2.1. Kevlar: 

Kevlar Aramid is an aromatic organic compound of carbon, hydrogen, oxygen, and 
nitrogen, and is produced by spinning long-chained polyamide polymers. With its material 
properties of low density, and high tensile strength, and low-cost, Kevlar produces tough 
and impact-resistant structure, however, the compressive properties of Kevlar laminates 
are low, due to poor coupling of resin matrices to the aramid fibers If compared to 
other organic fibers, the aromatic polyamide fibers of Kevlar, are characterized by their 
high tensile strength and high modulus 

There are two different types of Kevlar— Kevlar 29 and Kevlar 49 Out of the 
myriad of numbers of fiber materials that exist today, we chose a Kevlar type 49 weave for 
the majority of the wing construction primarily for its high strength and low weight It 
was also chosen for its failure quality Kevlar does not bond directly to the resin matrix 
This allows for the fiber to deform, yet still hold and maintain some structural integrity, in 
the case of an isolated impact Other fibers will bond directly to the resin, causing the 
overall material to be brittle In the same scenario these materials, if they fail, would fail 
catastrophically and shatter, leaving the aircraft structurally compromised Moreover, 
Kevlar 49 aramid fiber is preferred for the Banshee '94 structural design because it is 
characterized by low' density, high tensile strength and modulus. It is used in high- 
performance composite applications where light weight, high strength and stiffness, 
vibration damping, and resistance to damage, fatigue, and stress rupture are the main 
properties If Kevlar 49 is applied for reinforcement, the reinforced composites can save 
up to 40% of the weight of glass fiber composites at the same stiffness 

The stress-strain behavior of Kevlar 49 is linear to ultimate failure in tension at 
340 kips/in : and 1 8% elongation Usable strength of Kevlar 49 reinforced epoxy is 
about 4 times of 7075-T6 aluminum at less than half the density In addition, Kevlar 49 


61 



on a pound-for-pound basis is five times as strong as steel. The Figures 6.1a & 6.1b 
illustrates tensile stress/strain curves for tensile loading * * 



Figure 6.1a Specific tensile modulus. 10" irr 



Figure 6.1b Stress and strain curves of various fibers/ 
epoxy. 


( Reproduced from Airframe Structural Analysis, by Michael C.Y. Niu) 



Kevlar 49 is nonconductive and exhibits no galvanic reaction in contact with metal 
In addition it can be combined with glass or graphite in hybrid composites to take 
advantage of the unique properties of each fiber 

Epoxy resins are used to reinforce Kevlar Their high adhesion makes them useful 
in many high-performance applications in aircraft industry Epoxy resins have also 
excellent water resistance and low shrinkage during cure or vacuum bagging ( about 3%) 
Their conductive properties offer good fiber-matrix adhesion 

However, Kevlar does not cut well after it has been treated with a resin matrix 
When cut, it has a tendency to fray, and the edge will look somewhat like a shag carpet 
This would make it hard to cut the trailing edges of the control surfaces, which is done 
once the surface have been completed For these reasons, we have also decided to use 
carbon fiber in designated areas of the craft 

6.2.2. Carbon Fiber: 

The carbon fiber was also applied to areas where high strength was needed Since 
the area around the components in the upper surface have holes in them for access, it was 
decided that a layer of carbon fiber would be used around the holes to regain some of the 
strength lost The duct was also made of carbon since it would have no foam for support 
and rigidity 

To further the strength in the duct, joint, and component areas, two carbon fiber 
spars were inserted through these areas These extend to both rear landing gear 
placements to add support for their mountings 
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6.2.3. PVC Foam Core: 

The composite skin was built on a foam core. The purpose of this core is to act as 
a cushioning medium, to withstand the majority of the shear forces, and add to the rigidity 
of the wing 

There are two types of foam commonly used in composite structures Polystyrene 
foams and polyvinyl chloride (PVC) foams For the "Banshee 94", we opted to use the 
PVC foams for many reasons First, PVC is much stronger than polystyrene The foam 
has excellent tension, compression, and shear strengths that vary with the chosen density 
of foam This reduces the load that the skin is required to support, and allows for fewer 
skin layers increasing strength Another reason is its compatibility with composite skins 
It has better adhesion qualities with resins for a superior bond and stress communication 
PVC foams do not melt with resins, thus allowing good control of the wing shape Due to 
the PVCs' closed cell construction, the resin cannot permeate the foam, causing undesired 
weight, and thus maintaining a uniform density Lastly, PVC is easy to work with It can 
be hot wired and sanded to shape without crumbling. For the Banshee '94, we chose the 
Klegecell, type 45, 3 lb PVC We felt that this foam would be the lightest while 
maintaining the wings strength and rigidity with only one skin layer 

6.3 Methods of analysis: 

Today, there are many finite element methods which permit the efficient 
application of stress analysis : 

a) Abaqus and Ansvs are used for non-linear dynamics. While good for 

anisotropic materials, but they are poor for composites 

b) Argus is an excellent method for analyzing composites, non-linear post 

buckling and vax 
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c) I-Deas is a graphical method of analysis 

d) Nisa is good for doing an analysis on composite and vax 

e) Msc/Natran is excellent for structural analysis, but poor for some non-linear 
structures * . 

The structure groups chose the graphical method I-DEAS and the structural 
method of analysis, Msc/Nastran 

6.3.1. I-DEAS CODE: 

I-Deas is a finite element modeling code produced by Structural Dynamics 
Research Corporation To run the program, one needs a CAD drawing of the object to 
be analyzed After this the object has to be broken down to the desired number of nodes, 
or elements The properties of the materials to be analyzed are then entered the code in 
the form of elastic, and shear modulus The third step is to enter the forces and restraints 
to be applied to the model The results are displayed graphically with different colors 
representing different stress, strain, or deflection ranges 

For our wing, we decided to break the wing into four elements in the chord-wise 
direction, and twenty three elements along the span Since I-DEAS is not programmed to 
handle the complexities of composite structures, we were forced to make a number of 
assumptions First, since I-DEAS models a shell with an arbitrary thickness, or as a solid, 
we chose the model our wing as a shell The maximum thickness of the wing lies on the 
borderline between the definition of a solid and shell Since our main structural 
component, the Kevlar skin, is a thin sheet on either side of the wing, we felt that the shell 
with arbitrary thickness would give the best results. This would also cause the program to 
over estimate the stress, allowing for a worst case scenario 
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Since the model can be made of only one material for use of I-DEAS, our 
composite had to be converted to a single "materiar to enter the material properties 
Since the calculation for the ratio for Young's Modulus of the foam and Kevlas proved to 
be small it was felt that the elastic modulus of the PVC could be neglected However, for 
the shear modulus, the values for the PVC were entered since it was believed that it would 
be supporting the majority of these stresses 

Finally, the forces for level flight cruise conditions were applied to the model The 
vertical loading was the resultant forces from the elliptical lift distribution, component 
weight, the weight of the wing, and the winglet at the tip combined A graph of this load 
is shown in Figure 6 2 The drag loading was given by the Aerodynamic group, and is 
shown in Figure 6.3 The moment of 10 ft-lbs., applied at the quarter chord, was given by 
the Control group Finally, the wing root was restrained in all three axis and restrained 
from rotation A diagram of all forces applied to the wing can be seen in Figure 6 4 


fefcsKT 3 LfcVixlO ' 3 




Figure 6.2: Lift-Weight Loading 


Figure 6.3: Drag Loading 
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After running the program, the results were printed. The readings for XY plane. 
Von Mises, and maximum principal stresses can be seen in Figures 6.5 and 6.6. The 
maximum stress occurs in the Von Mises diagram which, in short, is the total of all 
stresses occurring at a given point This maximum value of 1 195 KSI falls well within the 
238 KSI of the Kevlar skin The diagram of the deflection, exaggerated in Figure 6.7 , 
shows a maximum deflection as being almost two thousandths of an inch. We felt that 
these values are acceptable for confirming the soundness of the structure 



Figure 6,4: Load Forces Applied to Wing in I-Deas 
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Figure 6.6: Von-Misis and Maximum Principle Stresses 
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Von-Mises Stresses Maximum Principal stresses 



Figure 6.7: Wing Deflection 
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6.3.2 Msc/Nastran Code: 

Msc/Nastran is a new analysis technique which generates a data file for one of the 
four types of problem s^ 

a) Plane stress analysis (Using CQUAD4 OR CQUAD8 elements). 

b) Beam (Using bar elements). 

c) Plate (Using CQUAD4 elements). 

d) Solid (Using CHEXA elements) 

Msc/Nastran uses a Finite Element Structural model, wherein the distributed 
physical properties of a structure are represented by a finite number of structure elements, 
which are interconnected at a finite number of grid points Loads are applied to grid points 
and for which displacement are calculated With the specific assumptions, the procedures 
for defining and loading a structural model for "Banshee'94",are generated through eight 
steps (see Appendix D1 ) 

1 Executive control deck- problem identification 

2 Case control deck - title ( Winfem Nastran Output, page!) 

3 Grid point definition - Grid point identifications and coordinates ( pages 3 to 6) 

4 Grid point constraints - Grid point identifications and constraint codes (p 3 to 
6) 

5 Elements - Elements identifications, property identifications ,and grid point 
connectives (pages 3 to 6) 

6 Properties - Material identifications and thickness (or section properties , if 
required) (pages 3 to 6) 

7. Material - Young Modulus, Poison ratio (pages 3 to 6) 
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8 Load grid point identifications and Cartesian components of forces (p. 3to 6) 

In comparing the results from Msc/Nastran to those from the I-Deas code, we find 
that the differences are great The Maximum Principal Stress for Msc/Nastran and I-Deas 
respectively, at the point of maximum value for I-Deas, is 3.983 lb/in 2 and 564 33 lb/in 2 . 
For the shear in the X-Y they are 0 0436 lb/in 2 and 113.62 Ib/in 2 . Finally, for the Von- 
Mises Stresses, we have 2 766 lb/in 2 and 1195 Ib/in 2 . The large differences in these 
values suggest that the different assumptions have a great impact on the results of these 
codes From these numbers we can also find that the I-Deas values may be higher than 
they should be To further investigate this, we performed an experimental test with a beam 
section 

6.3.3 Experimental Analysis: 

To check whether the assumptions used for the FEM codes were correct, it was 
decided to compare the strains and deflection of a composite test specimen with the data 
obtained of a computer simulation of the same specimen The specimen was cantilever 
beam 20 inches long, 2 inches thick, and 1 inch wide It was made of foam core with a 
layer of Kevlar laminated to the top and bottom surfaces. At the fixed end of this beam 
where the maximum strain would occur, a strain gauge was affixed Weight between zero 
to fifteen pounds were then hung from the free end, and the strain was measured with 1 5% 
experimental error (Fig 6 8) The deflection was also measured, and found to be near zero 
until the specimen failed The test beam withstood approximately 13.5 lbs before failing 
Upon failure, the lower surface Kevlar buckled, and the beam suffered massive 
deformation The beam was then simulated with the FEM codes and assumptions used for 
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the wing All forces were applied, including the distributed weight of the beam itself The 
results are shown in Figure 6 9 From this, we can conclude that the I-Deas codes are a 



Figure 6.8: Test Beam Dimensions Figure 6.9 Comparison of Test 


and I-Deas Data 


6.3.4 Comparison of analytical ethods with Msc/Nastran 
6.3.4. 1. VV eight/Strength 

Minimizing the weight of the aircraft has been a major consideration in the present 
design This led to the use of composite materials such as PVC, Kevlar, and Graphite, as 
well as new types of structures such as a geometrical block wing containing all loads 

From the des.gn point of view, with a given loading systems and the limitations on 
the outside dimens, ons (geometry) of the structure (if such limitations exist), a designer 
can figure out the optimum configuration of the internal structure. The loading system 
may include all the requirements as to breaking strength, deflection, fatigue etc .... 
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In determining the optimum design, the most important loading criterion is the 
ultimate strength Stiffness may also be required to modify the optimum structure and in 
some cases, may be very important as to outweigh strength completely. Because of the 
importance of providing adequate strength at minimum weight, it has become customary 
to evaluate materials and structures via their strength/weight ratios 

A strength/weight factor can be obtained by dividing the ultimate tensile strength 
of material by its density - 

Strength/Weight = o u + p = (a * F S) p (equation 6 1 a ) 
StrengthAVeight = (p * F.S) + ( A * p) (equation 6 1b) 

where 

o u Ultimate tensile strength 
p Densit\ 

F S Factor of Safety 
p Normal Force 
A Area 


Such a factor has the dimension of length It represents the length of vertical 
member, of a uniform section, at which the material fail in tension under its ow n weights 
Moreover, the factor is considered as the limit of maximum displacement at which the 
materials can withstand (Appendix D2) 

Factor ( Strength/Weight ) py^ = ] 52 inches 
Factor (Strength'Weight) PYC^Kevlar = ' 60 inches 
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6.3. 4. 2. Check Shear Stress and Bending Moments 

An additional method for checking shear stress and bending moments along 
the spanwise and the chordwise is based on the total area of the given geometrical tapered 
wing The total area of the whole wing(565 .57 in' ) is considered as a total load, which is 
assumed acting on the half wing in the Z direction and through the aerodynamic center of 
the airfoil section, and which is distributed along spanwise (see Appendix D3) 

6.3.4.3. Section Properties of the Composite Materials 

PVC foam "Q45" is considered as a homogeneous material with its Young 
Modulus(E), being a constant value However, the flexural rigidity which is defined as a 
product of a Young Modulus and the Moment of Inertia of each section (El), varies along 

the beam of varying its thickness In the case the wing is treated as a beam the product 

El takes an important role in analysis of deflection due to bending and torsion, therefore, 
the flexural rigidity must be calculated and their results are 

1 For PYC foam only 

EI(E329) = 922480 lbin : E 1 ( E32 7 ) = 330174 lbin : 

E1 ( e 328) = 512193 lbin ; E1( E326) = 150195 lbin : 

E1(E325) = 42279 Ibin' 

The procedure of flexural rigidity calculation is described and shown in the 
Appendix D4 

2 For the composites of PVC foam and Kevlar 

E1(E329) = 4460000 lbin : EI( E327 ) = 2050000 lbin : 

E1(E328) = 3240000 lbin' EI(E326) = 1180000 lbin' 

E1(E325) = 360000 lbin" 
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The procedure to calculate sections properties of the composites is described in (Appendix 
D4) 

With respect to aeroelastic effects, a tapered swept wing of a control surface is 
specifically important In the case, the deflection is due to bending, and torsion contributes 
considerably to the twist angle of the wing section The deflection at the tip on the elastic 
curve of the wing, treated as a beam, from a tangent to the elastic curve at the root is 
defined as the moment of the area under the M/EI diagram between points taken about the 
tipA 

S KT = J^(A/y 1 El)dy (equation 6.2) 

where 


S RT the spanwise deflection 

y the distance between centroid of elastic load (M/EI) and the tip 



\ 


Figure 6.10 Wing twist angle caused by swept. 

(Reproduced from Airframe Structural Analysis, by Michael N Y. Niu) 


76 



The result for the deflection due the bending moment of the PVC foam is 1.14 
inches, whereas, the deflection due to the bending moment of the composites of PVC and 
Kevlar decreases to the value of 0 0249 inches As its results is compared to the deflection 
of the composites obtained from the Nastran program (0 0284 inches), the value from 
hand-calculation is short off about 0 0035 inches This means that, as the commercial 

program is run, it produces the total deflection, being defines as 

$ + (equation 6 3) 

However, for the engineering calculation, the deflection due to the shear force is 
considered so small that its value can be neglected This is the reason why the value from 
the hand calculating is smaller than that obtained from Nastran program (see Appendix 
D5) 

Taking into account the limit of the factor of Strength/Weight of the materials for 
the composites ( I 60 inches), as well as comparing the results obtained from the Nastran 
program (maximum displacement - 0 0284 inches, including PVC foam and Kevlar) and 
from hand-ca!cu!ating(maximum displacement = 0 0249 inches), the structures group can 
conclude that the wing in the static case can withstand under its basic loads 

The concept of structural elastic axis is also applied for Torsion ( or chordwise 

deflection) The angular deflection due to torsion about the elastic axis between point A 

and B is given by the equation' ' 

B a - 6 B = j ( T / GJ )dy (equation 6 4) 

Where GJ is known as the torsional stiffness 


77 



Figure 6. 1 1 Wing chordwise deflection. 


(Reproduced from Airframe Structural Analysis, by Michael N Y. Niu) 

The values of torsional stiffness are obtained for each section as 

JG(E329) = 637200 Ibin' JG(£327)= 227610 lbin' 

JG(£ 328 ) = 334800 lbin' JG(£326) = 101520 lbin' 

JG(£325)= 26783 lbin : 

The above equations show that changes in airload due to structural 
deflection are due to torsional loads and that bending loads contribute more and more as 
sweep increases 

6.4. Construction Methods: 

Construction of the wing was begun by cutting the foam core This was done by 
first cutting the block of foam to the shape of the wing planform. Obtaining the desired 
airfoil shape was made easy by designing the wing in such a way that the midspan airfoils 
were linear extrapolations of the root and tip foils This meant that by making templates 
for the root and tip airfoils, affixing them to their respective ends of the foam block, and 
dragging the hot w'ire to coincide at specific chordlengths at each template, the midspan 
airfoils were automatically produced 
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When the core had been cut, the control surfaces were cut out Since the 
components were to build into the wings, small cavities were hollowed out for them as 
well as the channels for the controls rods and electrical wiring On the bottom of the wing 
cavities were hollowed out for the landing gear fixtures Into these were fitted blocks 
made of alternating layers of 1/32 inch plywood and 1/4 inch balsa The effect of these 
blocks are to dampen the shock of landing from the landing gear 

After fitting the control rods and electrical wiring, the foam received a thin layer of 
micro balloons on the surfaces that would be covered by the skin layer Micro-balloons 
are small, hollow spheres made of glass or phenolic resin that is mixed with the resin to 
reduced its weight while maintaining its strength The effect of this layer was to reduce 
weight of the aircraft by filling the small nooks on the surface of the foam with something 
that is lighter than untreated resin 

On top of this the Kevlar skin was added in a wet lay-up. and then vacuum- 
bagged Vacuum-bagging is the process in which the component is placed between layers 
of absorbitent material and then in an air tight bag Then a pressure of 25 inches Hg is 
applied to the component This process eliminates air bubbles in the resin that could cause 
stress build up in the skin as well as removing excess resin The duct was produced by 
molding layers of carbon fibers on a polystyrene form, and then melting out the 
polystyrene 

The wing were then connected to the duct by means of two carbon fiber rods that 
passed through the duct At the joints, we also applied gussets of Kevlar, 2 to 3 inches 
wide, for reinforcement We, then, applied another layer of micro-balloonmixture to fill 
the small "dimplies" that are caused by the weave in the skin. After this was completed, 
the components were mounted on balsa blocks in their respective compartments 
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Lastly, we mounted the engine and fan unit This was done by simply connecting the fan 
unit directly to the duct with eight screws, and rubber washers to dampen vibration 
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7. Conclusions and Recommendations 


Many difficulties were encountered during the design and construction of the Banshee 
’ 94 . In retrospect, many of these difficulties could have been either avoided or lessened, 
however these revelations were not obvious to the project team. 

As was expected, the series of calculations for the aerodynamic analysis of the Banshee 
’94 proved to be tedious and time-consuming. The use of a previously published vortex panel 
computer code proved to be more of a hinderauce than a help due to debugging and difficulties 
in deciphering the output. If such a code is to be used in the future, it would be best to know 
exactly how the code works and have a better than working knowledge of the language in which 
it is written. Another problem experienced by the Aerodynamics group is the wind tunnel 
testing apparatus. WPI does not have a working force-balance capable of measuring lift, drag, 
and moments on an airfoil simultaneously. An apparatus does exist, but it is not in working 
order. If the desire for such an apparatus is apparent, it is advised that work begin on making 
this apparatus operational early on in the project’s development. The Aerodynamics team 
concluded that the Banshee ’94 will fly at a cruise condition of angle of attack of 4 degrees and 
a cruise velocity of 33 miles per hour. 

Stability w-as one of the primary problems addressed in the design of the Banshee ’94. 
The group decided that the best way to handle the stability problems inherent to a flying wing 
was to oversize the control surfaces. In conjunction with the Aerodynamics group, a series of 
airfoils were selected which had a positive moment coefficient. This is essential in counteracting 
the instability of the flying wing. As experienced by the Aerodynamics group, computer codes 
proved to be a hinderance due to their difficulty in interpretation. Another problem addressed 
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by the Stability/Controls group was the implementation of the Electro-static Stability System. 
The existing literature available on the electrostatic stability system was rather minimal to the 
extent that several ambitious design implementations had to be undertaken. It is recommended 
that before an undertaking such as this, a better than working knowledge of electronics is 
acquired. From the analysis performed, the Banshee ’94 will maintain stable flight and will be 
able to maneuver without any major difficulty. 

The propulsion system that was used in the Banshee ’94 proved to have sufficient power. 
This helped insuring that enough power was available for take-off and maneuvering. It is 
recommended that future projects use the ducted fan engine assembly as their primary propulsion 
system. Testing of the engine proved to be a major problem due to an inappropriate force 
balance. Future projects should restrain the thrust gage in the wind tunnel before engaging the 
test, in order to ensure that no injuries or damages to the apparatus occur. 

The Structures group used two different finite element analysis codes to estimate the 
shear forces and displacements of the wing. Such codes are not entirely accurate due to the 
anisotropic nature of composites. The assumptions that had to be made when entering the data 
into the code proved to be difficult. The materials should prove to be more than adequate in 
handling the aerodynamic forces as well as the shocks from landing. 

The project team was very satisfied with the outcome of this project. The group found 
that the areas mentioned above were difficulties that could have been avoided had the team had 
the foresight to compensate for them. It is the view of this team that these areas be addressed 
early when tackling such a project. Further research and development may have provided a 
better design. More experimental data would have proven useful in determining some of the 
necessary values needed for calculations. 
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APPENDIX A2 


Lotus Spreadsheet for Drag and Lift Calculations 


c329 = 

C32 7= 

c325 = 

cg= 

Stot4wng 

b 

AR 

0.4148 
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0.2455 

13 . 6925 

0.73 

2.2098 

6.689337 

Stot* 

q= 

b= 
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Cmcg 
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329 
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325 


a. c. 
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13 
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-0.0072 

0.044 

-0.00457 

329 

327 
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Cd: 
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45 
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50 

305.9575 

0.987159 

0.226290 

0.306834 

0.021018 

24.585 
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STABILITY ANALYSIS 


B.l CONTROLS/STABILITY SPREADSHEET 
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B.2 INVISCID VORTEX PANEL CODE 


-13 


1 

f 

I 


PROGRAM VORCODE 


INTEGER Z,ZZ 
PARAMETER (Z=36) 
PARAMETER (ZZ=37) 


x ®( zz ) #YB(ZZ) ,X(Z) , Y(Z) , S(Z) , SINE(Z) , COSINE (Z) 
,CP ( Z > ' GAMA ( Z ) ,RHS(Z) ,CN1(Z,Z) ,CN2(Z,Z) 
DIMENSION CT2 (Z, Z) ,AN(ZZ, ZZ) , AT(Z, ZZ) , THETA ( Z ) ,CT1 (Z, Z) 


OPEN (UNIT= 11 , FILE='juju', STATUS= ' NEW ' ) 
OPEN (UNIT=12 , FILE='e327', STATUS= ' OLD ' ) 


M = 36 
MP1 = M+l 


DO 13 1=1 , MP1 

READ (12,*) XB (I) , YB (I) 

PI = 4.0 * ATAN (1.0) 
ALPHA = (4.) * PI/180 


DO 1 1=1, M 
IP1= 1+1 


X(I)=0.5* (XB (I) +XB (IP1) ) 

Y(I)=0.5 * (YB (I) + YB(IPl)) 

S(I)-SQRT( (XB(IPl) -XB(I) ) **2+ (YB (IP1) -YB (I) )**2) 

SIPS™ f (YB(IPl) -YB (I) ) , (XB(IPl)-XB(I) ) ) 
SINE (I) =SIN (THETA (I) ) 

COSINE (I) =COS (THETA (I) ) 


RHS (I) =SIN ( THETA ( I ) -ALPHA ) 
DO 3 1= 1 ,M 


DO 3 J= 1 , M 
IF (I.EQ.J) GO TO 2 

A=~(X(I) -XB(J) ) *COSINE(J) - (Y (I) -YB ( J) ) *SINE(J) 
B= (X (I) -XB ( J) )**2+(Y(I) -YB ( J) )**2 J 

C=SIN ( THETA ( I ) -THETA (J) ) 

D=COS (THETA (I) -THETA (J) ) 

E=(X(I)-XB(J) ) *SINE (J)-(Y(I) -YB ( J) ) *COSINE(J) 
F=ALOG (1 . 0+S ( J) * (S ( J) +2 . *A) /B) 

G=ATAN2 (E*S ( J) , B+A*S ( J) ) 

P=(X(I)-XB(J) ) *SIN (THETA (I) -2. *THETA(J) ) 

+ + ( Y (I) - YB ( J) ) *COS (THETA (I) -2 . *THETA ( J) ) 

Q-(X(I)-XB(J) ) *COS (THETA (I) -2 . *THETA ( J) ) 

+ “ (Y (I) -YB ( J) ) *SIN (THETA ( I ) -2 . *THETA ( J) ) 

CN2 (I , J) =D+. 5*Q*F/S ( J) - (A*C+D*E) *G/S(J) 

CN1 (I, J)=. 5*D*F+C*G-CN2 (I, J) 

CT2 (I, J)=C+.5*P*F/S( J) +(A*D-C*E) *G/S(J) 

CT1 (I, J)=. 5*C*F-D*G-CT2 (I, J) 


GO TO 3 


2 


CN1 (I , J) =-l . 0 


VD oo 


CN2 (I, J) *1. 0 
CT1(I, J)=0.5*PI 
CT2 (I , J) =0 . 5*PI 

3 CONTINUE 

DO 4 1=1, M 

i AN (1,1) =CN 1(1,1) 

AN (I ,MP1) =CN2 (I ,M) 

AT (I, 1) =CT1 (1,1) 

AT (I ,MP1 ) =CT2 ( I ,M) 

DO 4 J=2,M 

AN (I , J) =CN1 ( I , J) +CN2 (I , J-l) 

_4 AT (I , J) =CT1 (I , J) +CT2 ( I , J-l) 

AN (MP1 , 1 ) =1 . 0 
AN (MP1 ,MP1) =1.0 

DO 5 J=2,M 

_5 AN (MP1 , J) =0 . 0 

RHS (MP1) =0 . 0 


CALL CRAMER (AN, RHS, GAMA, MP1) 

DO 8 1=1, M 

V ( I ) =COS ( THETA ( I ) -ALPHA ) 

DO 7 J=1 , MP1 

V(I) =V (I) +AT(I , J) *GAMA ( J) 

CP(I)=1.0-V(I) **2 

WRITE (11,9) I , X ( I ) ,Y(I) , THETA ( I ) , S ( I ) , GAMA (I) ,V(I) ,CP(I) 
FORMAT ( 10X , 12 , F8 . 4 , F9 . 4 , F10 . 4 , F8 . 4 , 2F9 . 4 , F10 . 4 ) 

WRITE (11,10) MP1 , GAMA (MP1) 

10 FORMAT (10X, 12 , 35X, F9 . 4) 

END FILE (UNIT=11 ) 

END FILE (UNIT=12 ) 

CLOSE (UNIT=11) 

CLOSE (UNIT=12 ) 

STOP 

END 


SUBROUTINE CRAMER (C, A, X,N) 

INTEGER ZZ 
PARAMETER (ZZ=37) 

DIMENSION C(ZZ,ZZ) ,CC(ZZ,ZZ) ,A(ZZ) ,X(ZZ) 

DENOM=DETERM ( C , N ) 

DO 3 K=1 , N 
DO 1 1=1, N 
DO 1 J=1 , N 
CC (I , J) =C ( I , J) 

DO 2 1=1, N 
CC ( I , K ) =A ( I ) 



X (K) =DETERM ( CC , N) /DENOM 

RETURN 

END 


FUNCTION DETERM (ARRAY, N) 

INTEGER ZZ 
PARAMETER (ZZ=3 7) 

DIMENSION ARRAY (ZZ,ZZ) , A (ZZ , ZZ) 
DO 1 I-1,N 
DO 1 J=1 , N 

A (I , J) =ARRAY (I , J) 


M=1 

K=M+1 


DO 3 I=K,N 

RATIO=A ( I , M) /A (M, M) 

DO 3 J=K, N 

A (I, J)=A(I, J) -A(M, J) * RATIO 

IF (M.EQ.N-1) GO TO 4 

M=M+1 
GO TO 2 
DETERM=1 

DO 5 L=1 , N 

DETERM=DETERM* A ( L , L) 

RETURN 

END 
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Data from the Climb Rate vs. Velocity 
and calculations for Acceleration 


Velocity 

(mph) 


RC 

(ft/s)( 1 00%) 


0 


32.16 


39.02 


45.3 


51.15 


55.39 


60.08 


62.55 


65.67 


66.14 


67.06 


64.87 


63.18 













RC 
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13.169251 
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18.85035 
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0 
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0 
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Accel (75%) 

Accel (50%) 


(100%) 



0 

37.56675 

28.17506 

18.78337 

0.2423 

35.30657 

25.91488 

16.5232 

0.3786 
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APPENDIX D2: STRENGTH / WEIGHT 

FACTOR 



Using the equations 6 1 and 6 2 in the chapter 6.3 (section 6.3.3 1), PVC factor 
(Strength/Weight) is calculated 


Factor (Strength/Weight) = o u ! p 


Area of the wing = 565 57 in" 

Strip load = 1 lb 

F S = 15 unmanned aircreaft 

p (density of PVC) = 3 Ib/ft 5 = 0 0017361 1 Ib/in' 

Then 

Factor for PVC = (lib* I 5)/(565 57in : * 0 00173611 Ib/in ' ) 
Factor foi PVC = 1 52 in 

The Factor of composite PVC and Kevlar is obtained as the following 
Density of PVC = 0 001 7366 Ib/in ' 

Density of Kevlar = 0 052 lb/in' 

Ratio of Density = p K v , / p 

Ratio of Density = 0 052 /0.00173611=29 952 
Therefore 

Factor (Strength/Weight) = o u I p d „ c „„, r 

Factor (Strength/Weight) = ( 1 5* 1 045 1 6 )/( 565.57 * 0.0017361 1 ) 
Factor (Strength/Weight) = J 60 in 



APPENDIX D3 : CALCULATING 

WING SHEARS AND MOMENT FOR 
ONE UNIT LOAD CONDITION 



Calculation of the wing shear V z and the wing moments M x and My due to the 
total unit , distributed on a half wing load of 565.57 lbs, acting upward in Z direction 
and applied at aerodynamic centers of the wing sections goes through the following 
steps 

Step 1 The distance from the wing root section, called station y (Column I ) 

Step 2 The wing chord length of each section (Column 2) 

Step 3 Ratio of Lift coefficient is assumed to be unity (Column 3) 

Step 4 The w ing running load per inchof span at each station point The area of a half 

wing (565.57 square inches) was treated as a total unit load, then, the ainning 
load per inch at any station equals the wing chord length at that station 
(Column4) 

Step 5 Average running load (Column 5) 

Step 6 The distance between stations (Column 6) 

Step 7 Strip load ( Column 7) 

Step 8 The location of the strip load which is at the centroid of a trapezoidal 

distributed load whose end values are given in column 4 The centroid locations 
are determine by Table A3 4 of Chapter 3 (Column 8) 

Step 0 Shears (Column 9) 

Step 10 The product of shear and delta y (Column 10) 

Step 1 I The product of strip load and arm to the centroid (Column 1 1 ) 

Step 12 The value of moemnt (Column 12) 

Step 13 The distance from aerodynamic center to the y reference axis (Colun 13) 

Step 14 The strip load distance from y axis (Column 14) 

Step 1 5 Bending moment for each strip load on the y axis (Column 1 5) 

Step 16 The bending moment at the various stations equals to the summation of the trip 

moments as one progresses from station 43 .5 to zero. 


(Ref 6? 
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APPENDIX D4 : 

SECTION PROPERTIES 
OF 

THE COMPOSITE MATERIALS 



The product EJ of each wing section (E329, E327, E326, E326, E325) is calculated 
through the following steps: 


Step I 
Step 2 
Step 3 

Step 4 

Step 5 

Step 6 
Step 7 

Step 8 

Step 9 

Step 10 

Step 1 I 


The geometry of each section is known. 

The reference points of coordinates system is chosen 

The cross section area is divided into available partitions and numbered from 1 
to 1 2 or 1 3 ( which depend on each cross section area) . ( Column 1 ) 

The dimensions of each partition are measured in z and y directions, called 
delt z and deit y (Columns 2 and 3) 

The area of each partition is calculated, and after that areas of all partitions are 
summed up (Column 4) 

The centroid Z c „ and Y ctJ for each partition is calculated (Columns 5 and 6) 
The product of area and the distance from the centroid in z and y directions are 
found, and after that they are summed up (Columns 7 and 8) 

The product of area and square of Z ca and square Y C g are calculated , and 
summed up ( Columns 9 and 10) 

The Inertia at the centroid of each partition are found in the z and y direction 
and summed up ( Columns 1 1 and 12) 

The product of area and Z C g and Y C g are calculated, and summed up 
(Column 13) 

Torsional constant is defined as 


- For a rectangular section 



- For a quarter of an ellipse section: 

7,, = (U16)KZ x(K : +Z 2 ) 

- For a rectangular section 


J„ = (K 4 V3)/80 

where 

K = Vy : + Z 2 

(Column 14) 

(Reference Engineering Mechanics Statics, by R C Hibbeler. section Centroid and 
Center Gravity, page 457) 

Step 12 Z(BAR) and Y(BAR) are the centroid components of an object in Yand z 
directions, and determined by the centroid equations: 


Z = 
Y = 




(Reference Engineering Mechanics Statics, by R C Hibbeler, section Centroid and Center 
Gravity, page 420) 

Step 12 IzzOyv and I zv are calculated by the equations 


i z/ = £az : + £i ,-z : a 
4 -y : a 

l vz = AYZ-AYZ 

Step 14 The product of flexural rigidity is defined as 

EI U = Z pVGir.-. * 4 

Step 15 The product of torsional stiffness is 





Step 16 The torsional stiffness applied at the centroid is given by the equation 


./ = ./ - AY Z 

t, o 

The whole procedure is shown in the tables 1,2. 3,4, and 5 for each section E325, E326, 
E327, E328 and E329 

(Reference Introduction to Aerospace Structural Analysis, by David h. Allen and Walter 
E. Haisler. section of Determination of Modulus Weighted Section Properties 4 3 2 3, 
page 1 77) 

Finally, the graph represents the flexural rigidity and the torsional stiffness versus the 
stations of the wing 
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APPENDIX D5 : THE RESULTS 

FOR CHECKING DEFLECTION FOR 

PCV CORE 
AND 

FOR COMPOSITES-PVC & KEVLAR 



The procedures for calculating the vertical static deflection for the PVC wing core are 
going through the eight steps 

Step 1 The wing station for each airfoil section is determined. (Column 1) 

Step 2 The change in distance along the spanwise (Column 2) 

Step 3 The bending moments are taken from the Column 16, Table of Calculation of 
wing shear and wing moments. Appendix C (Column 3) 

Step 4 The flexural rigidity is imported from Tables E325, E326, E327, E328, E329, 
Appendix D (Column 4) 

Step 5 The ratios of the bending moments and the product of flexural rigidity are 
calculated (Column 5) 

Step 6 The area of each section is calculated (Column 6) 

Step 7 The deflection of each section is determined. (Column 7) 

Step 8 The total is deflection for the whole wing is the sum of the deflections from the 
root to the tip 

Similarly, the deflection for the wing, made of PVC and Kevlar, is calculated After 
that, a graph shows the deflection versus the wing stations 
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